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Abstract 


This report describes the development of an improved long-wave radiometer (pyrgeometer) 
for deployment on ships and buoys. Standard pyrgeometers use a thermopile to measure the 
temperature gradient between the receiver surface and the instrument case, and thus infer 
the receiver temperature and incident radiation. The key design change employed in the 
new radiometer is to remove the thermopile and replace it with a small, glass-encapsulated 
thermistor to measure the receiver temperature directly. To prove the concept, a prototype 
radiometer was built and calibrated. It was then deployed outside for a period of a week on 
the roof of the Clark Laboratory (Quissett Campus, Woods Hole Oceanographic Institution) 
to demonstrate the feasibility of the new concept. Data from the prototype were compared 
to those from a pair of standard radiometers. The intercomparison shows that the prototype 
performed surprisingly well. It was able to capture all the variability observed by the stan¬ 
dards with only a small bias. The next step in the design process, which has been funded 
by the National Science Foundation, is to build a rugged version of the prototype that can 
be deployed in the field. 

1 Introduction 

Recent climate research programs are demanding an increased accuracy of air-sea flux es¬ 
timates. The surface long-wave (infrared) radiation fluxes have not commonly been mea¬ 
sured at sea, and improvement of the long-wave radiation instrumentation is needed to meet 
these demands. Improved meteorological measurements and aerodynamic bulk formulae 
have brought net errors in short-wave radiation and latent and sensible heat fluxes down to 
approximately 10 W m~ 2 . Conservative error estimates for typical long-wave radiation flux 
measurements are 10 W m -2 at night and 40 W m -2 during the day. Thus the long-wave 
radiation is the source of the largest uncertainty in air-sea surface energy budgets. 

Long-wave radiation is a difficult measurement to make from a ship, since the ship itself 
is often a source of infrared radiation and internal heating effects in the radiometers make the 
instruments nearly useless unless care is taken to shade the instrument from direct sunlight. 
The Eppley Precision Infrared Radiometer (PIR) has been the standard in the atmospheric 
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research community for over 20 years and is often deployed on land and aircraft. Surprisingly, 
few users fully understand the instrument or are aware of (or just accept) its limitations. 
And even fewer have access to calibration facilities for these instruments. The Woods Hole 
Oceanographic Institution now has one of only a handful of calibration facilities in the world 
(Payne and Anderson, 1999). The expertise gained in setting up this facility has placed us 
in a unique position to begin design improvements. 

In this report, we describe an effort to develop an improved long-wave radiometer 
(pyrgeometer), for deployment on ships and buoys, using our improved understanding of the 
instrument and modern measurement techniques. The new design overcomes some of the 
limitations of the currently available instrumentation, and will lead to a reduction of the 
errors in measuring surface long-wave radiation. 

In our redesign of the pyrgeometer, we endeavored to develop an instrument that met 
the following criteria: 

• is easily calibrated and the calibrations remain stable even under unfavorable condi¬ 
tions; 

• has few components and is easily fabricated; 

• is not significantly affected by short-wave radiative heating effects of the radiometer 
body; and 

• is rugged and reliable enough to withstand a 12-month deployment on a surface buoy. 

Standard radiometers, like the Eppley PIR, employ a thermopile to measure the re¬ 
ceiver surface temperature. This technique was required to increase the instrument’s sensi¬ 
tivity. Eliminating the thermopile is a revolutionary idea, and its good sense is the major 
finding presented here. The thermopile has been in use for decades and the research com¬ 
munity has been in a bit of rut. We found no need to use a thermopile in the instrument 
since the required sensitivity could be achieved with inexpensive glass-encapsulated thermis¬ 
tors. Removing the thermopile will reduce costs of fabrication and improve the measurement 
capability. We envision the community embracing this technique in the coming years. 

A review of theory of pyrgeometer operation is presented in Section 2. In Section 3, 
we discuss our considerations for redesigning the Eppley PIR. Our experiences in designing 
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and fabricating a prototype and validating the new measurement technique are in Section 4. 
We successfully demonstrate that the prototype not only works but performs remarkably 
well. The Appendix includes some information on using empirical formulae for estimating 
long-wave radiation at the ocean surface. 


2 Background 

The net long-wave flux at the sea surface is the difference between the downward (LW 4-) 
and upward radiation (LW t) • 

LWU = LWt -LWl (1) 

LW tl is a function of sea surface skin temperature T s *, sea surface emissivity, e, and the 
profiles of temperature and humidity in the lower atmosphere. The upward flux is given by 

LW t = eaT^ (2) 

where a is the Stefan-Boltzmann constant. Empirical formulae assume that in a cloud- 
free regime, the net long-wave radiation can be described as a function of the sea surface 
temperature, T s , the air temperature, T a , and the near-surface water vapor pressure, e a . 
These formulae are effective at predicting the net long-wave radiation in clear sky. The 
difficulty is in parameterizing the effect of cloud cover using observable quantities. Cloud 
cover observations are taken by shipboard observers who are trained to estimate cloud base 
height and cover by eye. An alternative approach is to use measured incoming solar insolation 
and compare it with the theoretical clear sky values to estimate the cloud cover. This method, 
however, does not yield any information about cloud height. (See Fung et al., 1984, for a 
review of empirical formula.) 

The preferred method of estimating net long-wave flux is a “hybrid measurement sys¬ 
tem” where the incoming long-wave radiation at the sea surface is measured and the outgoing 
long-wave radiation is estimated from sea surface temperature measurements. The formula¬ 
tion is 

LWU(0) = LWt (0) - LWL (0) 

= e ss cjTj - (1 - e ss )LWl (0) (3) 
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where e ss is the sea surface emissivity (0.985 ± 0.010), T ss is the sea surface temperature, 
and LWi (0) is the measured incoming long-wave radiation flux at the sea surface. 

2.1 Pyrgeometer theory 

In the ideal pyrgeometer, the downwelling radiation is estimated from the balance of absorbed 
radiation to emitted radiation from a surface. The heat budget for the surface is simply 

H — R ne t 

= aRl-e 0 aT s 4 (4) 


where H is the heat gain to the surface from conduction or convection, R net is the net 
radiation on the surface, T s is the receiver temperature, Rl is the incoming radiation from 
the environment, a is the absorption coefficient of the surface (a = e 0 = 1 for a black body), 
and e 0 is the surface emittance. In the ideal case, H is zero, T s is measured, then the Rl is 
calculated. In practice, no instrument can be built with this simple heat budget. 

The first difficulty is that the receiving surface cannot be supported and T s cannot be 
measured without creating a path for conductive heat loss. Thus there must be a case that 
supports the receiving surface and dome. If the heat loss to the case is only conduction, it 
can be written as 

k(T s - T c ) (5) 

where A; is a heat transfer coefficient for the conductive heat loss between the surface and 
the case. (Note that this term should be as small as possible. External heating of the case 
may lead to thermal gradients along the path between T s and T c which are not just a simple 
function of the conduction of heat away from the receiving surface. This would lead to a 
poor estimate of the conductive heat loss from the surface to the case.) Adding in this term, 
the heat budget then becomes 

k(T s - T c ) = aRl -eaT s 4 (6) 


Note that k will depend on materials and geometry of the case. 
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The second difficulty is designing a filter that restricts transmitted radiation to long¬ 
wave only. This is usually implemented with a hemispherical silicon dome with an interfer¬ 
ence filter coating on the inside. The dome transmittance must be a function of wavelength. 
It will be assumed here that it has one value in the short-wave band and another in the 
long-wave band. The transmitted radiation is 

RX = r$ w SWl + t% w LWI (7) 

where SW X and LW X are incident short-wave and long-wave radiation fluxes and rf w and 
r d w are the transmittance in each band. The dome unfortunately also acts a source and 
reflector of radiation. Thus the net radiation of the receiver surface is 

7?net ~ Rl + 7?T (8) 

where 

i?T = Co <jT s 4 + (1 — eo)-R4- (9) 

and 

Rl= T $ w SW l +t£ w LW l +pRt +eaT d A (10) 

where p and e are the long-wave reflectance and emittance of the dome and Td is the dome 
temperature. Note that in the above equation, it is assumed that there are no short-wave 
emissions from either the dome or the receiving surface. Taking note of the fact that (1— p) — 
e + r, we can solve for LW X directly 

LWi = oT A + E—4 Ir^ HTs - T c ) + —£w(T s 4 - T/) - (11) 

T d e o r d t^ 

(I) (II) (III) (IV) 

The four terms are (I) the black body emission from the surface, (II) the heat lost due to 
conduction to the case, (III) the radiative heat flux from the dome, and (IV) the short-wave 
leakage through the dome. This is the same as Equation 7 of Fairall et al. (1998) except for 
the inclusion of the short-wave transmittance. The LW X could theoretically be estimated 
from measured SW !• T s , T c , and T d . 
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2.2 The Eppley PIR 


The standard pyrgeometer, designed about 30 years ago, is the Eppley Precision Infrared 
Radiometer (PIR). This instrument employs a silicon hemisphere dome which provides a 
nominal passband of approximately 4-50 /an, isolating terrestrial long-wave radiation from 
the solar short-wave radiation. The radiation is sensed using a non-wavelength selective 
thermopile detector located beneath the dome and a thermistor located in the case. The 
thermopile consists of 60 thermocouples in series to increase the sensitivity. The thermopile 
measures a temperature difference rather than an absolute temperature. The hot junction 
of the thermopile forms the receiving surface, and the cold junction is in thermal contact 
with the case. The voltage output of the thermopile is 

so ^out = (T s ~ T c ) = AT (12) 

where s 0 is the thermopile sensitivity. One of the thermistors, T c , is placed as close as 
possible to the thermopile cold junction. Another thermistor, T d , is attached to the inside of 
the dome. Substituting these temperatures into the above equation, the incoming long-wave 
is 

LWl = a(T e + AT)* + kAT + [(T. + AT)* - T,‘] -%SW{ (13) 

T d e 0 T d T d 

(I) (II) (HI) (IV) 

To make this equation look like the formulation by Albrecht and Cox (1977) 

L = E(ci + C 2 r c 3 ) + €q<t T c 4 - ka(T d 4 - T c 4 ) (14) 

the following assumptions must be made: 

r! W = 0 

AT T c so that (AT + T c ) 4 « T c 4 + 4T C 3 AT (15) 


Then the downwelling long-wave equation reduces to 


LWl = °T C A + s-'V,out 


t NM! + 4 / lt 

T d w e o V 


£d \ 

LW I 


crT c 3 


T; 


+ 




'*4 

c 


T?) ( 16 ) 
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[This formulation is different from that of Alados-Arboledas et al. (1988, Eq. 2). They have 
a term T c 3 term which is not multiplied by AT. It is not clear where this term could come 
from since it must be derived from the Taylor series expansion.] 

In the standard Eppley pyrgeometer configuration, the contribution of the last term in 
Eq. 16 is ignored. The Eppley sensor uses a compensation circuit to combine the thermopile 
output voltage, the case thermistor resistance, and a reference battery voltage. The voltage 
out of the unit is directly proportional to LW | (see Fig. 1). This configuration was used 
in both the vector-averaging wind recorders (VAWR) and improved meteorological (IMET) 
PIRs deployed on the WHOI surface mooring during TOGA COARE (Tropical Ocean Global 
Atmosphere Coupled Ocean-Atmosphere Response Experiment). When the units are cali¬ 
brated, there is a potentiometer that should be adjusted to tune the compensating circuit 
for zero error at typical deployment temperatures. We have not tested this but it is unlikely 
that this was done. (For full description of the errors from the compensation circuit, see 
Olivieri, 1991.) 

The current version of the IMET PIR measures the three outputs, thermopile voltage, 
dome and case thermistor resistances, separately, and computes temperatures and total long¬ 
wave irradiance by Eq. 16. Also, the IMET PIR has an all-aluminum case with no radiation 
shield whereas the stock Eppley PIR is constructed with brass, stainless steel and aluminum 
with a radiation shield. See Payne and Anderson (1999) for details. At the time of the initial 
PIR design, Eppley had been building pyranometers based on a clever use of a thermopile 
for some time. It probably was not feasible, given the thermistor technology of the time, 
to design a pyrgeometer which used thermistors exclusively. Use of the thermopile was a 
reasonable choice. 

2.3 Calibration coefficients 

Equation (11) can be simply rewritten using a coefficient for each of terms II, III and IV 
such that 


LWi = crT s 4 + A(T S — T c ) + B(T s a — T d 4 ) — C SW l (17) 

(I) (II) (III) (IV) 
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Figure 1: Eppley Precision Infrared Pyrgeometer circuit schematic. 







Table 1 

Suggested values for Eppley PIR radiative terms 




Dickey et al. 
(1994) 

Olivieri 

(1991) 

long-wave transmissivity of the dome 

Td 

0.27 ±0.01 

0.34 

short-wave transmissivity of the dome 

~sw 

T d 

0.0098 ± 0.005 

0.003 

long-wave reflectivity of the dome 

Pd 

0.4 

— 

short-wave reflectivity of the dome 

P s d w 

0.99 

— 

long-wave emissivity of the dome 

td 

0.3 ±0.1 

— 

short-wave emissivity of the dome 

c hW jsw 

6 5 € d 

0.0 

— 

long-wave emissivity of the receiving surface 

6 

0.91 ±0.01 

0.985 

reflectivity of the receiving surface 

p LW^ p SW 

0.02 ±0.01 

— 


where A, B and C are coefficients determined by calibration. The coefficients A and B can 
be determined directly with a black body cavity using the procedure described by Payne and 
Anderson (1999) who do not evaluate C. The B term is directly related to the Albrecht and 
Cox (1977) k term. The A coefficient, however, is not directly relatable to the C\ and c<i terms 
for Albrecht and Cox (1977, their Eq. 15). Note that, because of the difficulties of applying 
the interference coating to the dome interior, the long-wave transmittance and emissivity 
and the short-wave of the dome vary from one dome to another leading to a variation in 
values of B. 

Dickey et al. (1994) suggest values of all the radiative constants (See Table 1). They 
employ the PIR equation written as 

LWi = 0.91 crT; 4 + s; 1 l/ out (0.976R + 8.08aT c * 3 ) ± 1.11(T/ - T c 4 ) ± 0.03675IT; (18) 
(I) (II) (III) (IV) 

The value of B is determined from calibration, and Dickey et al. do not attempt to propose 
a value for it. Note that the coefficient of 0.91 is compared with 1.0 multiplying term (I) in 
Eq. 17. Using this coefficient, a value of T c = 300°K corresponds to a difference in LW 4- of 
41 W m~ 2 . 
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Dickey et al. arrive at short-wave transmittance by assuming that most of the so¬ 
lar effect seen by Alados-Arboledas et al. (1988) is from solar radiation being transmitted 
through the dome. Olivieri (1991) proposed a much smaller short-wave transmittance. He 
assumes the transmission of the silicon dome is negligible below 2.5/xm. Since solar irradiance 
is negligible above 5/un, one only needs to consider a small band of short-wave radiation. 
Olivieri’s estimates would lead to noon time errors in LW i of only 1.2 W in -2 for a SW l 
of 400 W m -2 as opposed to 14.78 W m~ 2 from Dickey et al. Olivieri also suggests a much 
larger long-wave transmittance of r d = 0.34. 

Olivieri proposes that the major solar error signal is from overheating of the dome. 
This is consistent with Albrecht and Cox (1977) who present observations of the heating of 
the dome and its effect on measured downwelling long-wave radiation. They find that this 
error is nearly a direct function of observed short-wave radiation. They suggest an empirical 
correction for direct solar heating of the pyrgeometer dome using observed SW |. 

A SW I 

ka(T* 4 - Tf) = a SW i + (19) 

This formulation uses the time-lagged observed short-wave radiation to account for the ther¬ 
mal mass of the pyrgeometer. This may be important for the fast response observations 
they took from an aircraft, however, this term will be small in the moored observation time 
series. They suggest a solar heating coefficient, a, of 0.0311 (i b = 0.0666) which is nearly the 
same as Dickey et al. It is important to note that this coefficient is proportional to their k. 
It appears that the formulation of Dickey et al. is flawed. Their formulation should either 
have the dome temperature correction term or the empirical correction using the observed 
but not both. It appears that the short-wave transmittance is indeed small and can be 
largely ignored while the dome heating term must be accounted for in long-wave calculations. 

2.4 Thermal gradients and other sources of measurement error 

Thermal gradients within the pyrgeometer are another source of error. This is not accounted 
for in the theory nor given much attention. Foot (1986) has identified thermal gradients 
in the case and dome. Measurement errors will result since the case temperature may 
not be representative of the cold junction side of the thermopile. In addition, thermal 
gradients across the dome mean the dome thermistor does not accurately represent the 
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Table 2 

Summary of uncertainty for Eppley PIR 


long-wave measurements 




(I) 

term 

(II) 

(III) 

total (W m~ 2 ) 

T c 

±0.100°K 

±0.6 


±1.3 

±1.9 

AT 

±0.005°K 

±0.03 

±2.0 

±0.1 

±2.0 

T d 

±0.200°K 



±2.6 

±2.6 

A 

4% 


±3.1 


±3.1 

B 

4% 



±0.8 

±0.8 

Total uncertainty 

±0.6 

±3.7 

±3.0 

±5.0 


average temperature of the dome (Philipona et al., 1995). There have been attempts to 
measure the dome temperature at three different locations around the dome simultaneously, 
but there are no documented attempts at measuring the case temperature gradients directly. 

Another source of error is the thermistors themselves. Payne and Anderson (1999) 
report that the PIR thermistors must be calibrated since some depart significantly from 
the manufacturer’s specification, particularly the dome thermistors. This is not typically 
done however. There is no technique available for calibrating the thermopile sensitivity 
directly, and its small output voltages place a very low noise and high accuracy demand on 
the amplifier and A/D sampling electronics. There also remains some uncertainty in the 
calibration coefficients. 

A su mm ary of the uncertainties in the Eppley PIR measurements is given in Table 2. 
These assume that T c = 300°K, T d = 301°K, AT = — 0.2°K, A = 390, and B = 2.8. It 
appears that for times longer than a week, the PIR accuracy may approach ±6 W m~ 2 as 
validated by intercomparisons (Fairall et al., 1998; Payne and Anderson, 1999). However, 
the instantaneous accuracy can be much worse, close to ±30 W m -2 . 


15 



3 Rethinking the pyrgeometer design 

The Eppley PIR uses a thermopile to measure the receiver surface temperature. This tech¬ 
nique was required to increase the instrument sensitivity, however, it makes the instrument 
susceptible to thermal gradients across the case and there is significant conductive heat loss 
from the receiving surface to the case through the thermopile. In addition, the size of the hot 
junction accounts for only 17% of the area below the dome, which results in large amounts of 
emitted and reflected radiation from the case being absorbed by the detector. Finally, the in¬ 
strument cases are made of brass and stainless steel which have poor thermal characteristics 
and are susceptible to corrosion in the marine environment. 

Here, we reexamine the PIR and seek to improve on the design in three specific ways. 
First, new glass-encapsulated thermistors and associated electronics allow for the direct 
measurement of the detector temperature with the required sensitivity and would replace 
the thermopile. This allows for better isolation of the detector from the case and for easier 
calibration. Secondly, the case would be constructed to minimize the conductive heat loss 
of the receiver to the case and reduce temperature gradients in the case. Thirdly, the size of 
the detector would be enlarged to nearly 100% of the available area under the dome which 
would increase the detector sensitivity. 

3.1 Replacing the thermopile 

The thermopile has been a tool in radiation measurement for decades. These devices are 
expensive to fabricate, and an absolute calibration is usually not available. But is the pre¬ 
cision of temperature differences the thermopile provides really needed in the pyrgeometer? 
Can we eliminate the thermopile altogether and measure the temperature of the receiver 
surface directly using a glass-encapsulated thermistor and still achieve measurement uncer¬ 
tainties of 5 W m -2 or better? This would allow for a cheaper-to-build and easier-to-calibrate 
pyrgeometer. 

With a glass-encapsulated thermistor, we conservatively can measure temperatures 
to ±0.010°K and often much better. Using Eq. 17, a 0.010°K error in T s translates into a 
(I) 0.1 W m” 2 error in the black body emittance term, (II) 3.9 W m“ 2 error in the conductive 
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heat loss term, and (III) 0.1 W m~ 2 error in the dome heating term (assuming A = 390 and 
B = 2.8). The only term in the pyrgeometer equation that requires the high sensitivity is 
the conductive heat loss term. This result was not intuitive to us. 

If the pyrgeometer can be redesigned such as to reduce the conductive heat loss 
term (III), the sensitivity achieved with the thermopile is not needed. This would require 
increasing the dome transmissivity in the long-wave band, increasing the emissivity (eo) of 
the receiver surface, or decreasing the conductive heat loss coefficient, k, to the case. The 
emissivity of Eppley-Parsons black lacquer, used on the receiver of the PIR, is 0.985, so 
there is not much room for improvement. The transmissivity of the dome in the long-wave is 
approximately 0.35. It is unlikely that the dome long-wave transmissivity could be increased 
without increasing the short-wave transmissivity, which would cause other problems. If it 
could be increased, it would also reduce the dome reflectance and emittance terms and re¬ 
duce the dome heating errors. However, we did not explore this option since it is beyond 
our expertise at this time. The only domes available to us are those made by Eppley. 

We can, however, attempt to reduce the conduction of heat from the receiver to the 
case. In fact, doing away with the thermopile gives us the flexibility to do this. In the 
thermopile, the thermocouple wire is wound around an aluminum T-shaped support (see 
Fig. 2), and the cold junction of the thermopile must be in good thermal contact with 
the case through the arms of the T. This assures a good conductive heat path. Without 
the thermopile, we are free to insulate the receiver from the case, effectively reducing the 
coefficient k. What we found in our subsequent testing is that this term can, indeed, be 
readily reduced by nearly a factor of 4. In our prototype, discussed in the next section, the 
coefficient A is 110 while for most Eppley PIRs this number is near 400. 

We note here, and again later, that although reducing the heat loss to the case is 
important, just as important is still being able to accurately estimate term (II) in the in¬ 
strument since there will always be some conduction to the case. We found that using case 
materials with low heat conductivity supported large temperature gradients in the case and 
that no meaningful T c measurement could be made to effectively quantify term (II). 

The reduction of the need for high sensitivity is not the only advantage to replacing the 
thermopile. By reducing the magnitude of the heat loss term (III), we reduce the importance 
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Figure 2: Upper panel: Photograph of thermopile from Eppley PIR. Lower panel: Photo¬ 
graph of Eppley PIRs, left, stock Eppley PIR; right, IMET PIR. 
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of the calibration and place more of the contribution to the long-wave measurement from 
term (I), which has no calibration coefficient. 

Finally, removing the thermopile removed the necessity of high gain amplification as 
well. This circuit has proved susceptible to R/F interference. Using the thermistors requires 
just a simple resistance bridge, which is easy to design, is stable, and requires low power. 
The associated electronics will be much simpler and robust. 

3.2 The dome heating correction 

Reconsider the dome heating, term (III). In the Albrecht and Cox (1977) formulation, their 
k equals the e d /r d w or B in the equations. They report that k may be as large as 4.0 for a 
stock Eppley sensor. Olivieri (1991) also suggests that this coefficient is about 4.0. Due to 
the restraint that e d + r d + p d = 1, it is not possible to reach a coefficient of 4.0 using either 
of the suggested values by Dickey et al. (1994) or Olivieri for long-wave emittance. 

In the above formulation it has been assumed that the area of the thermopile surface 
is equal to the total area under the dome. This is not the case in the Eppley pyrgeometer 
(see Fig. 3). The dome is attached to a removable cover. In the IMET version, the cover 
is aluminum with an exterior coating of white paint. The VAWR version has a cover made 
of stainless steel. In both configurations, a lip of this material, A 2 , is exposed to the area 
under the dome. The case in both sensors is made of aluminum and also has a surface, Ai , 
that is exposed under the dome. The true thermopile surface, A s , is a fraction of the total 
area under the dome. (The geometry under the dome also includes a thermister, which is 
glued to the inside of the dome. The area of the glue holding the dome and thermistor is 
Aq-) An accurate assessment of the heat budget on the thermopile surface would include 
the emittance and reflectance of all these surfaces. 

For example, consider the B(T C 4 — T d A ) term (III). Taking the long-wave emittance 
from the additional surfaces, ei, 62 , and tc, and assume the surface temperatures on A\, A 2 , 
and Aq are equal to T d . Using a direct and reflected path for the emitted radiation, the 
coefficient becomes 


(1 + PdXuAi + 62^2 + egAg) 
Td T d (Ai + A2 + A s ) 


( 20 ) 



3.0 



dimensions in cm 


c) schematic of area under dome 



Figure 3: Schematic of stock Eppley PIR. 
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This formulation allows for larger possible values of B. Taking the emissivities as 


^stainless 

^aluminum 

£glue 


0.6 (slightly oxidized) 

0.1 

0.9 


the normalized areas 

^Totai — + A 2 + A s ; (X\ = j4i/Axotai = 0.242 ; 02 = 0.6 ; clq = 0.1 

and the values suggested by Dickey et al. lead to the estimates 

B = 3.6 for stock Eppley (VAWR) and B = 2.0 for IMET . 

Thus, the different materials used in the case will lead to large differences in the sensitivity 
to the dome temperature. The IMET PIR is made of aluminum and should be less sensitive 
to temperature differences between the dome and the case than the stock Eppley PIR. 
The constant B has been estimated by Dickey et al. for four different IMET sensors. The 
individual values range from about 1.4 to 2.2. Fairall et al. (1998) show a table of values of 
B from all of our calibrations of IMET PIRs. The range of values of 3.2 to 4.4 is consistent 
with the above calculations. 

A simple experiment was conducted to test the hypothesis that the cover material and 
geometry has a significant effect on the heat budget estimation on the thermopile surface. An 
Eppley PIR [S/N 27185; 4.19 mv/(W m" 2 )], an IMET PIR [S/N 27926; 4.19 mv/(W m" 2 )], 
and an Eppley 8-48 short-wave sensor were placed side by side in a shaded room. Both PIRs 
were equipped with the internal circuit and no dome temperature measurements were made. 
The long-wave sensors were allowed to reach equilibrium and yielded the same estimates of 
LWl- The shade was removed, allowing a direct path for the solar insolation for 140 seconds, 
then the shade was replaced. Two typical runs are shown in Figure 4. Eppley specifications 
give a response time of 2 seconds. The response observed has two time scales. The first 
is the rapid response of the thermopile reacting to the additional long-wave radiation. The 
response is short (4-8 seconds) because the mass of the thermopile is small. After the initial 
rapid rise, the output continues to rise for the duration of the solar insolation. The VAWR 
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Figure 4: Example of radiative heating in stock IMET and VAWR type PIRs. In Run 8, an IMET and a VAWR 
PIR, along with a pyrgeometer, were exposed to sunlight for 140 sec then covered. Note rapid initial response of the 
long-wave as the thermopile responds, then the gradual increase as the dome and cover begin to heat up. The VAWR 
has a much larger error at the end of 140 sec. In Run 11, the covers from the VAWR and IMET were switched. The 
larger error is found to track with the VAWR cover, not with the case. 





output, however, rises much faster, climbing to 7 W m -2 above the IMET output at the 
end of the solar insolation. When the solar insolation is removed, the outputs drop rapidly 
again and quickly reach a value that is larger than the initial output by the amount equal 
to the increase from the slow response over the time of insolation. This same pattern was 
observed even if the sensor positions were switched or various shields were placed over the 
sensor cases leaving only the domes exposed. 

The next step was to switch the covers and domes from the sensors so that the IMET 
had the stainless top from the VAWR and vice versa. Then the experiment was repeated. The 
same pattern was observed on various conditions except this time the output from the IMET 
rose above that of the VAWR. The excessive output due to insolation followed the cover, not 
the sensor. If there were a pin hole in one of the domes or the short-wave transmissivity were 
different between the domes, the initial response between the two instruments would be of 
different magnitudes. This is because the thermopile temperature will respond quickly to the 
increase in infrared radiation. There is a small difference (2-3 W m -2 ) between the two covers 
just after the initial rise in the output (10 seconds after the start of insolation). The longer 
time scale response must be related to heating of the cover and the dome that slowly increases 
the radiation incident on the thermopile surface. This is not a rigorous experiment and 
should be redone measuring the individual thermistors and thermopile outputs separately. 
However, it is clear that the dome heating term tracks with the PIR dome and top and not 
the thermopile and case. 

The effect of dome heating on the receiver can be reduced by increasing the area of 
the receiving surface under the dome. Like reducing the conductive heat loss term, reducing 
the dome heating term will lessen the dependence of measurements on calibration, putting 
more of the dependence on term (I), the black body radiance term, which is independent of 
calibration. Expanding the receiver area is easily addressed in the redesign that removes the 
thermopile. 

4 Demonstration of concept 

In this section we describe our efforts to build a prototype pyrgeometer based on the ideas 
described in the previous section. The prototype was built to be easily modified so we 
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could make changes in the lab to explore the thermal characteristics and responses. The 
instrument was not designed for field deployment but rather as a proof of concept. However, 
we did deploy the prototype alongside an IMET PIR on the roof of the Clark Laboratory 
for two weeks. From this deployment we do successfully demonstrate that the prototype not 
only works but performs remarkably well. This effort has proved the concept of removing 
the thermopile completely. The next step would be to develop the instrument for long term 
deployment. 


4.1 The prototype 

A drawing of the prototype designed and built is shown in Figure 5. The body is made 
of al uminum and is painted with a bright white paint. The receiver surface is aluminum, 
0.11 mm thick, painted with Parson’s optical black lacquer. Its diameter is 2.00 cm as 
compared with the diameter of the dome base, 2.67 cm. Attached with thermal conducting 
epoxy to the middle of the underside is a 0.36 mm-diameter thermistor. The receiver disk is 
cemented to the tops of three 2-56 threaded rods, located near its perimeter, which project 
8.4 mm above an aluminum plug insert. This insert is screwed into the inside of the body. 
The stainless steel rods provide a convenient low thermal conductance mounting for the 
receiver disk. A thermistor is potted into the body and another is potted into the receiver 
mounting plug insert at the base of one of the receiver support rods. The silicon dome was 
purchased from Eppley and is identical to the domes used in the Eppley PIR. 

4.2 Calibration 

The prototype pyrgeometer was calibrated using the same technique and equipment we use 
for the PIR (Payne and Anderson, 1999). Six runs are made, each one lasting 12 minutes. 
In three of them, data are recorded for 1 minute while the PIR is in equilibrium at room 
temperature, for 1 minute while the PIR is suspended upside down over water at 50°C, and 
for 10 minutes while the PIR is mounted in a cold black body cavity at, successively, 0.1, 5, 
and 10°C. In the other three runs, data are recorded for 1 minute at room temperature, for 
1 minute while it is suspended over ice, and for 10 minutes while the PIR is mounted in a 
warm black body cavity at, successively, 30, 40, and 50°C. There is at least 1 hour between 
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Figure 5: Drawing of pyrgeometer prototype: (b) insert. 
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Figure 5: Drawing of pyrgeometer prototype: (c) overall sketch. 
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runs to allow the PIR to return to room temperature equilibrium. One of the warm-to-cold 
runs is shown in Fig. 6. Note the sudden rise of the plate temperature and the slower rise 
of the body temperature followed by a reversal when the pyrgeometer is put in the black 
body cavity. Where the curves cross the second time, T s = T c and term II in Eq. 17 is zero. 
Equation 17 then reduces to 

LWi= aT s + B(T S 4 - T/) (21) 

In later versions of the pyrgeometer when we measure dome temperature, we will determine 
the B constant. Since LW i is known from the measured temperature of the black body 
cavity, this will yield six independent determinations of B. Knowing B, we will compute A 
using the full Eq. 17 and the last data record from each of the six runs. In the prototype, 
without the measurement of the dome temperature, we compute six values of A from the 
last record neglecting term III. 

4.3 Field inter comparison 

Figures 7-10 show data from an 11-day comparison of the prototype pyrgeometer with two 
Eppley IMET-style PIRs on the roof of the Clark building, a site with a clear view of the 
whole sky. Short-wave radiation from an Eppley PSP is overplotted to show where that is 
significant. All sensor outputs were recorded on a Campbell CR7 data logger. 

The three pyrgeometers track very well except during daytime periods of high short¬ 
wave radiation when the lack of a dome correction in the prototype causes its values to be 
anomalously high. Except in these high short-wave radiation periods, the prototype agrees 
with the Eppley PIRs to 10 W m -2 . 

Figure 8 shows the components and total long-wave flux as computed from Eq. 17 
without the short-wave correction term. Figure 9 shows the difference between the prototype 
and one of the Eppley PIRs. Except for three very short periods, the differences are of order 
10 W m~ 2 or less. Figure 10 is a plot of prototype-PIR difference vs total PIR long-wave 
flux. The large excursions are principally during periods of high short-wave flux. There is a 
small dependence of the difference on the magnitude of the flux. 
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and body temperatures, respectively. 
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Figure 7: Total flux from two Eppley PIRs and the prototype pyrgeometer in an 11-day field test vs time. Overlay 
of short-wave flux. 
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Figure 10: Plot of difference between prototype and one Eppley PIR total long-wave flux for field test vs PIR total 




We found it quite encouraging that the comparison showed such small differences for 
a prototype intended to prove the concept. 


5 Conclusions 

Our study has shown that it is possible to build a pyrgeometer to measure downwelling 
long-wave radiation without the use of a thermopile. The new instrument is simpler to 
construct and calibrate, and has a simpler theory of operation, than standard thermopile- 
based pyrgeometers. The comparison of data from the prototype and standard instruments 
demonstrates that the concept clearly works. 

We have completed the first step in the development of a new meteorological instru¬ 
ment. The prototype pyrgeometer was able to demonstrate the concept but there still re¬ 
mains a lot of work before we have a field-deployable instrument capable of reliably collecting 
scientific-quality data in the field. 

The next step is to construct a field-deployable verison of the new pyrgeometer design. 
The prototype did not have a thermopile dome, and thus the dome heating errors could 
not be accounted for in the measurements. This must be included in the next version. In 
addition, the instr um ent, should be rugged enough to withstand the hazards of deployment 
in the marine environment and functionally similar to our existing inventory of standard 
pyrgeometers. Once such an instrument is constructed, it should be deployed next to a 
standard instrument during several field experiments to gain confidence in its reliability and 
accuracy. 


34 



Appendix 


Clark et al. formulation 

The net long-wave radiation formula of Clark et al. (1974) is commonly used and uses the 
same basic parameterization as competing formulations. (Note the formula often referred to 
as Clark et al. is based on work by Berliand and Berliand, 1952.) 

LWU= ™T S A [0.39 - 0.05(e a ) 1/2 ] F(C ) + 4 eaT s 3 (T s - T a ) (22) 

where the cloud cover factor is F(C) = (1 - bC 2 ). T s and T a in the formula are given 
in degrees Kelvin. The vapor pressure, e a (in mbars), is calculated using Buck’s (1981) 
formula for saturation vapor pressure of pure water corrected for salt water which used the 
measured relative humidity, barometric pressure and sea surface temperature. The cloud 
cover parameter 6 is a function of latitude: 

50° N 40°N 30°N 20°N 10°N 5°N Equator 

b 0.73 0.69 0.64 0.60 0.56 0.53 0.51 

The cloud cover, C, is the cloud cover fraction from 0.0 to 1.0 where 0.0 corresponds 
to a clear sky and 1.0 is overcast. The clear sky LW t! is larger than that of a cloud-covered 
sky with the same T s , T a , and e a . 


Cloud cover estimation 


The cloud cover can be estimated from the measuring incoming short-wave radiation. Using 
the cloud factor formula from Kimball (1928) 


C = 1.41 


SWobs I 

SWcsi 


(23) 


where 5W 0 bs l is the observed incoming solar radiation, S'Wcs i is the theoretical clear sky 
value for incoming solar radiation. The SWes I is calculated using the formulae given in the 
Smithsonian Meterological Tables (List, 1951): 
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5W CS ; = 


^WdifFuse + ^W^direct 

0.5^ (0.91 - a csc a ) sin a 
0.5^|(0.91 - a 080 a ) sin a 


+ 


Jo 


a csca sin a 


L r z 


( 24 ) 


where J 0 is the solar constant, a is the sun’s inclination, r is the radius vector of the earth, and 
a is the atmospheric transmission coefficient (ATC). The values of z and r are determined by- 
calculating the local solar time and using the formula found in “The Astronomical Almanac, 
1986,” US Naval Observatory, page C24. (Also see “Almanac for Computers 1988,” Nautical 
Almanac Office, US Naval Observatory.) The atmospheric transmission coefficient must be 
estimated using the SW ohs 4, during times of clear sky. 


a — 


2,W ob s|r- 2 

Jo sin a 


sin a 


-0.91 


(25) 


The ATC is a function of the sun’s zenith but is modeled here as a constant for all angles. 
Thus, values estimated at dawn and dusk will be different than those estimated at noon. 

There are many pitfalls in using the SW ohs l to estimate the cloud cover. When the 
sun is near the horizon, a cloud far away from the observation site may block the solar 
insolation yielding a high cloud-cover estimate even if there are no overhead clouds. The 
opposite may also occur when clear sky near the horizon leaves a direct path for. the solar 
insolation but the overhead sky is nearly completely cloud-covered. Thus the cloud-cover 
formulation of Kimball (1928) is not expected to yield accurate cloud-cover estimates on 
a minute-to- min ute (or even hour-to-hour) time scale. Another difficulty is that the cloud 
cover can only be estimated during daylight hours. The best that can be hoped for is that 
the estimated cloud cover averaged over the daylight hours will reflect a daily average cloud 
cover. 
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